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We  propose  a  series  of  icosahedral  matryoshka  clusters  of  A@B12@A20  (A  =  Sn,  Pb;  B  =  Mg,  Zn,  Cd),  which 
possess  large  HOMO-LUMO  gaps  (1.29  to  1.54  eV)  and  low  formation  energies  (0.06  to  0.21  eV/atom).  A 
global  minimum  search  using  a  genetic  algorithm  and  density  functional  theory  calculations  confirms  that 
such  onion-like  three-shell  structures  are  the  ground  states  for  these  A21B12  binary  clusters.  All  of  these 
icosahedral  matryoshka  clusters,  including  two  previously  found  ones,  i.e.,  [As@Ni12@As20]3_  and  [Sn@ 
Cu12@Sn20]  12~,  follow  the  1 08-electron  rule,  which  originates  from  the  high  /,,  symmetry  and  consequently 
the  splitting  of  superatom  orbitals  of  high  angular  momentum.  More  interestingly,  two  magnetic 
matryoshka  clusters,  i.e.,  Sn@Mn12@Sn20  and  Pb@Mn12@Pb20,  are  designed,  which  combine  a  large 
magnetic  moment  of  28  pB,  a  moderate  HOMO-LUMO  gap,  and  weak  inter-cluster  interaction  energy, 
making  them  ideal  building  blocks  in  novel  magnetic  materials  and  devices. 

The  first  evidence  for  the  formation  of  bare  anionic  clusters  of  the  post-transition  elements,  particularly  those 
of  groups  14  and  15,  was  obtained  by  Zintl  and  co-workers1-4  in  the  1930s  from  the  potentiometric  titrations 
of  post-transition  elements  with  alkali  metals  in  liquid  ammonia.  However,  definitive  structure  determina¬ 
tions  of  these  clusters  became  possible  only  in  the  1970s,  when  Kummer  and  Diehl  reported  the  first  structurally- 
authenticated  Zintl  anion  cluster  of  Na4Sn95.  The  original  post-transition  element  clusters  were  empty,  i.  e.,  they 
contained  no  interstitial  atoms  in  the  center  of  the  cluster  polyhedron.  However,  subsequent  synthetic  studies  led 
to  the  discovery  of  clusters  containing  interstitial  transition  metals.  Examples  of  structurally  characterized  such 
clusters  based  on  10-vertex  outer  polyhedra  include  the  anionic  indium  clusters  M@In1010‘  (M  =  Ni,  Pd,  Pt)6  in 
the  intermetallics  K10In10M,  the  anionic  lead  clusters  M@Pb102  in  [K(2,2,2-crypt)]2[M@Pb10]  (M  =  Ni,  Pd,  Pt)7, 
Fe@Sn103'  in  [K(2,2,2-crypt)]3[Fe@Sn10]8,  the  cationic  pentagonal  antiprismatic  bismuth  cluster  Pd@Bi104+  in 
Bi14PdBr16  (  =  [Pd@Bi10]  [BiBr4]4)9,  and  the  centered  pentagonal  prismatic  clusters  M@Ge103'  (M  =  Co10,  Fe11). 
Interstitial  metal  atoms  are  also  found  in  the  icosahedral  M@Pb122-  clusters  (M  =  Ni,  Pd,  Pt)7,12.  Furthermore, 
Al3+  has  been  encapsulated  in  Pb122-,  Pb102_  and  Sn102_  clusters,  which  have  been  characterized  by  mass 
spectroscopy  and  density  functional  theory  (DFT)  calculations13,14. 

Attempts  to  expand  the  knowledge  of  transition  metal  derivatives  of  post-transition  element  clusters  led  to  the 
discovery  of  even  more  complicated  structures.  Of  particular  interest  is  the  cluster  [As@Ni12@As20]3_,  isolated  as 
its  m-Bu4P+  salt  and  structurally  characterized  by  X-ray  diffraction15.  This  cluster  anion  has  a  perfect  icosahedral 
(Ih)  symmetrical  three-shell  matryoshka  doll  structure  consisting  of  an  outer  As20  regular  dodecahedron  encap¬ 
sulating  an  Ni12  icosahedron,  which  in  turn  encapsulates  an  As3-  trianion  (Figure  1).  Beyond  the  filled  d10  shells  of 
the  arsenic  and  nickel  atoms,  this  system  has  a  total  of  108  valence  electrons,  among  which  40  skeletal  electrons16 
correspond  to  a  “magic”  number  of  a  jellium  sphere17  with  lS2lPfTD102S2lF142P6  molecular  orbital  configuration. 

Initially  the  icosahedral  matryoshka  [As@Ni12@As20]3-  structure  with  108  valence  electrons  was  considered  as 
an  anomaly,  unique  with  an  unusual  combination  of  specific  electronic  and  steric  properties  of  its  arsenic  and 
nickel  components.  However,  in  2011  Stegmaier  and  Fassler18  reported  a  completely  analogous  [Sn@Cu12@ 
Sn20]12”  icosahedral  matryoshka  doll  anion  in  the  M12Cu12Sn21  intermetallics  (M  =  Na,  K).  This  [Sn@Cu12@ 
Sn20]12~  anion  is  exactly  isoelectronic  with  the  [As@Ni12@As20]3_  anion  with  108  valence  electrons  beyond  the 
filled  d10  shells  of  the  tin  and  copper  components.  Similar  bonding  models  appear  to  be  applicable  to  both  the 
[As@Ni12@As20]3~  and  [Sn@Cu12@Sn20]12-  systems.  Furthermore,  the  observation  of  the  icosahedral 
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Figure  1  j  Atomic  structure  of  an  icosahedral  matryoshka  cluster  of 
C@B|2@A2o  (left)  and  its  schematic  plot  (right). 

matryoshka  doll  structure  in  both  the  soluble  [As@Ni12@As20]3~ 
anion  and  the  [Sn@Cui2@Sn20]12-  anion  in  the  A12Cui2Sn21  inter  - 
metallics  suggests  that  such  structures  are  particularly  favorable  in 
clusters  where  108  total  valence  electrons  are  available  from  the 
elements  composing  the  three  layers  of  the  matryoshka  structure. 

Parallel  to  the  condensed  phased  clusters19  such  as  those  discussed 
above,  gas-phase  clusters  have  been  intensively  studied  during  the 
past  three  decades  since  they  exhibit  many  fascinating  physical  and 
chemical  properties  that  depend  on  their  size,  geometry,  and  com¬ 
position20.  Interestingly,  some  specific  clusters  with  appreciable 
stability  mimic  the  chemical  behavior  of  elemental  atoms  in  the 
periodic  table  and  thus  can  be  regarded  as  “superatoms”21-24.  In  these 
superatom  clusters,  the  electronic  states  are  delocalized  over  the 
entire  cluster  with  spatial  shapes  resembling  the  atomic  orbitals. 
The  corresponding  energy  eigenvalues  can  be  grouped  into 
atomic-like  shells  (IS,  IP,  ID,  2S,  IF,  2P,  ...),  depending  on  the 
degeneracy  and  spatial  symmetry  of  the  molecular  orbitals.  Unlike 
the  atoms,  the  physical  and  chemical  properties  of  superclusters  can 
be  tailored  through  selection  of  size  and  composition,  making  them 
very  promising  building  blocks  for  new  materials25,26. 

The  concept  of  “superatoms”  was  originally  proposed  by  Khanna 
and  Jena  in  their  pioneering  paper21.  For  example,  an  Al13“  cluster 
with  filled  IS2,  IP6,  ID10,  2S2,  IF14,  and  2P6  shells  exhibits  strong 
resistance  to  oxidation27.  Meanwhile,  with  an  electron  affinity 
(EA=3.57  eV)22  only  slightly  lower  than  that  of  the  Cl  atom 
(3.62  eV),  a  neutral  Al13  cluster  was  considered  as  a  superhalogen 
by  Bergeron  and  co-workers27,  though  this  nomenclature  is  not  con¬ 
sistent  with  the  original  definition  of  superhallogen  by  Bolydyrev  and 
Gutsev28,  which  defined  superhalogen  as  a  cluster/molecule  whose 
electron  affinity  is  larger  than  that  of  any  halogen  atom.  Moreover,  a 
neutral  Al14  cluster  behaves  like  an  alkaline  earth  atom29. 

Usually,  a  superatom  cluster  is  stable  because  it  has  filled  electron 
shells  with  paired  electrons,  corresponding  to  a  nonmagnetic  state. 
The  concept  of  magnetic  superatom  was  initially  introduced  by 
Kumar  and  Kawazoe30,  who  predicted  Mn@X12  clusters  (X=Ge, 
Sn)  as  icosahedral  magnetic  superatoms  with  a  high  magnetic 
moment  of  5  pB  and  a  large  highest  occupied  molecular  orbital- 
lowest  unoccupied  molecular  orbital  (HOMO-LUMO)  gap  of 
~1.1  eV.  By  appropriately  combining  the  localized  d  states  with 
magnetism  and  the  delocalized  superatom  orbitals  that  stabilize 
the  entire  cluster,  Khanna  and  co-workers  recently  designed  a  series 


of  magnetic  superatoms31,  such  as  VCs832,  VNas-33,  Mg8Fe34, 
Ca8Fe35,  MnCa936,  MnSr937,  ScK12  and  ScCs1238.  All  of  them  pos¬ 
sess  enhanced  stabilities,  substantial  FIOMO-LUMO  gaps  (usually 
about  0.4  ~  0.7  eV),  and  large  magnetic  moments  on  the  trans¬ 
ition  metal  atom  (up  to  5  pB). 

Stimulated  by  the  above-mentioned  progresses,  here  we  propose  a 
series  of  icosahedral  matryoshka  clusters,  i.e.,  A@B12@A20  (A  =  Sn, 
Pb;  B  =  Mg,  Zn,  Cd,  Mn).  Analogous  to  [As@Ni12@As20]3~  and 
[Sn@Cu12@Sn20]12-,  all  of  these  clusters  show  high  stability  and 
alignment  of  superatom  orbitals.  In  particular,  the  Sn@Mn12@Sn20 
and  Pb@Mn12@Pb20  clusters  possess  a  giant  magnetic  moment  of  28 
(J.B,  which  can  be  retained  in  cluster  assemblies.  Our  results  thus 
extend  the  scope  of  icosahedral  matryoshka  clusters  and  provide 
novel  building  blocks  for  cluster-based  materials  and  devices  exhib¬ 
iting  fascinating  novel  magnetic  and  optical  properties. 

Results  and  Discussion 

The  binding  energies,  formation  energies,  HOMO-LUMO  gaps,  and 
optical  absorption  gaps  for  icosahedral  matryoshka  A@B12@A20 
clusters  (A  =  Sn,  Pb;  B  =  Mg,  Zn,  Cd)  are  summarized  in  Table  1. 
For  a  given  cluster  of  AxBy,  its  formation  energy  is  defined  by:  Et 
(AxBy)  =  E(AxBy)  -  xE(A)  -  yE(B),  where  E(AxBy)  is  the  energy  of 
AxBy  cluster,  E(A)  or  E(B)  is  the  energy  per  atom  of  the  constituent 
element  A  or  B  in  the  standard  states,  x  or  y  is  the  number  of  A  or  B 
atoms  in  the  cluster.  Detailed  information  on  geometry  parameters 
and  on-site  charges  for  these  clusters  are  given  in  Table  SI  of 
Supplementary  Information. 

All  the  clusters  explored  here  show  reasonably  low  formation 
energy  (typically  between  0.1  and  0.2  eV/atom),  indicating  that  they 
are  easy  to  form  thermodynamically.  For  comparison,  we  consider 
several  38-atom  elementary  metal  clusters  with  a  fcc-like  truncated 
octahedron  structure39.  Higher  formation  energies  are  obtained,  i.e., 
0.225  eV/atom  for  Pb38,  0.281  eV/atom  for  Sn38,  0.576  eV/atom  for 
Mg38,  0.497  eV/atom  for  Zn38,  and  0.385  eV/atom  for  Cd38  clusters, 
respectively.  Using  global  optimization  combined  with  DFT  calcula¬ 
tions,  Ferrando’s  group40  has  recently  achieved  a  series  of  three-shell 
high-symmetry  matryoshka  clusters,  including  a  Ni13@Mg20  that  is 
isostructural  to  the  clusters  considered  here.  However,  the  formation 
energy  of  this  Ni13@Mg20  cluster  is  as  high  as  0.468  eV/atom  from 
our  DFT  calculation.  The  more  favorable  formation  energies  of  the 
present  icosahedral  matryoshka  superatoms  clearly  demonstrate  the 
stabilization  effect  of  electronic  structure,  which  will  be  discussed  in 
detail  later. 

Note  that  the  lengths  of  Pb-Pb  (3.317  to  3.578  A)  and  Sn-Sn  (3.176 
to  3.463  A)  bonds  listed  in  Table  SI  are  longer  than  the  equilibrium 
bond  lengths  of  the  hollow  Pb20  (3.073  A)  and  Sn20  (2.891  A)  cages 
from  our  DFT  optimization.  Therefore,  the  enlarged  Pb20  and  Sn20 
cages  are  stabilized  by  the  encapsulated  C@B12  icosahedron  (C  =  Sn, 
Pb;  B  =  Mg,  Zn,  Cd),  as  previously  found41  for  [As@Ni12@As20]3C  It 
is  also  interesting  to  compare  the  structural  parameters  of  Sn@Zn12@ 
Sn20  with  its  isoelectronic  counterpart,  i.e.,  [Sn@Cu12@Sn20]12-.  The 
theoretical  Sn-Sn  bond  lengths  in  Sn@Zn12@Sn20  of  3.278  A  are 
comparable  to  that  of  3.290  A  for  [Sn@Cu12@Sn20]12-  computed 
with  the  same  level  of  theory.  Note  that  the  experimental  Sn-Sn  bond 


Table  1  |  Binding  energies  (Eb),  formation  energy  (Ef),  HOMO-LUMO  gaps  (EHl),  optical  absorption 
transition  orbitals  for  the  icosahedral  matryoshka  clusters  A@B]  2@A2o  (A  =  Sn,  Pb;  B=Mg,  Zn,  Cd) 

gaps  (Eog) 

,  and  the  lowest  allowed 

Cluster 

Eb  (eV/atom) 

Ef  (eV/atom) 

EHL(eV) 

E„g  (eV) 

Lowest  allowed  transition  obitals 

Sn2iMg12 

2.79 

0.09 

1.38 

1.93 

HOMO- 

->  LUMO+1 

Sn2 i Zn ! 2 

3.43 

0.16 

1.47 

1.94 

HOMO- 

-»  LUMO+1 

Sn2iCdi2 

2.60 

0.21 

1.29 

1.84 

HOMO- 

-»  LUMO+1 

Pb21Mg12 

2.52 

0.06 

1.50 

1.74 

HOMO- 

2  — »  LUMO 

Pb2iZrn2 

2.97 

0.18 

1.54 

1.83 

HOMO- 

->  LUMO+1 

Pb21Cdi2 

2.28 

0.16 

1.27 

1.72 

HOMO- 

->  LUMO+1 
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lengths  of  [Sn@Cu12@Sn20]12-  in  the  condensed  phase  obtained  by 
X-ray  crystallography  range  between  3.076  to  3.133  A18. 

According  to  the  on-site  population  analysis  (see  Supplementary 
Table  SI),  the  B  atoms  in  the  intermediate  shell  donate  some  charge 
to  the  A  atoms  in  the  exterior  shell,  while  the  central  A  atom  gains 
more  charge,  ranging  from  0.245  to  1.362  electrons  (by  Mulliken 
definition).  The  trend  of  charge  transfer  remains  the  same  for  dif¬ 
ferent  systems,  suggesting  that  all  these  clusters  share  the  same  pat¬ 
tern  of  chemical  bonding.  The  direction  of  charge  transfer  can  be 
easily  understood  by  the  electronegativity  differences,  that  is,  the 
Pauling  electronegativities  of  Pb  (2.33)  and  Sn  (1.96)  are  higher  than 
those  of  Mg  (1.31),  Zn  (1.65),  Cd  (1.69),  and  Mn  (1.55). 

As  shown  in  Table  1,  the  HOMO-LUMO  gaps  of  the  A21B12 
clusters  with  B  =  Mg,  Zn,  Cd  range  from  1.27  to  1.54  eV,  comparable 
to  those  of  [As@Ni12@As20]3~  (1.44  eV)  and  [Sn@Cu12@Sn20]12- 
(1.39  eV)  using  the  same  computational  scheme.  Previous  DFT  cal¬ 
culations  yielded  similar  band  gaps  for  [As@Ni12@As20]3_ 
(1.44  eV41)  and  [Sn@Cu12@Sn20]12-  (1.34  eV18).  The  high  symmet¬ 
ries  of  these  clusters  brings  about  multiple  degeneracy  of  the  frontier 
molecular  orbitals  and  severe  selection  rules  in  optical  transitions. 
Our  TD-DFT  calculations  further  reveal  that  the  lowest  singlet 
excitation  (HOMO-to-LUMO)  is  optically  forbidden  due  to  the 
LUMO  and  HOMO  symmetries.  The  first  optically- allowed  excita¬ 
tion  (i.e.,  optical  absorption  gap  Eog)  and  the  associated  energy  levels 
are  summarized  in  Table  1.  The  computed  Eog  by  TD-DFT  ranges 
from  1.72  to  1.94  eV  (i.e.,  the  visible  light  region),  making  them 
potentially  useful  as  nanoscale  building  blocks  for  optoelectronic 
devices25,41. 

As  aforementioned,  the  experimentally  found  icosahedral 
matryoshka  clusters,  i.e.,  [As@Ni12@As20]3_  and  [Sn@Cu12@Sn20]12-, 
are  exacdy  isoelectronic  with  108  valence  electrons  beyond  the  filled 
d10  shells  (for  As,  Ni,  Sn,  Cu).  The  same  108-electron  rule  applies  to  all 
binary  A@B12@A20  clusters  considered  here,  noting  that  Mg  has  no  d 
electrons.  However,  108  does  not  explicitly  belong  to  any  existing  magic 
numbers  for  various  models,  e.g.,  the  shell  model  within  the  jellium 
approximation23,  the  Wade-Mingos  rules42,  the  spherical  aromatic 
model43,  the  octet  or  eighteen-electron  rule23.  Therefore,  we  explore 
the  origin  of  the  “magic  number”  of  108  for  these  spherical  multi-shell 
clusters  as  a  guide  to  design  other  matryoshka  clusters  exhibiting  high 
stabilities  and  sizeable  electronic  gaps. 

For  simplicity,  we  start  with  Sn@Mg12@Sn20  that  does  not  involve 
any  d  electrons  in  the  valence  shells  of  the  component  atoms.  The 
energy  levels  and  spatial  distribution  for  the  Kohn-Sham  molecular 
orbitals  of  Sn@Mg12@Sn20  are  presented  in  Figure  2  and  Figure  3, 
respectively.  According  to  the  energetic  sequence  and  the  nodal 
shape  of  the  orbitals,  we  can  identify  a  series  of  superatom  orbitals 
in  the  sequence  IS,  IP,  2S,  ID,  IF,  1G,  2P,  3S,  2D,  1H,  2F,  3P,  II 
corresponding  to  the  54  lowest-lying  valence  molecular  orbitals  (108 
valence  electrons)  up  to  and  including  HOMO.  Among  them,  IS,  IP, 
2S,  ID,  IF,  1G,  2P,  3S,  2D,  3P  orbitals  are  completely  filled.  In  the 
icosahedral  ligand  field44,  the  2F,  1H  and  II  superatom  orbitals  are 
split  into  two,  three  and  four  components,  respectively.  The  lowest  II 
orbital  (gg)  is  occupied  as  the  HOMO  of  superatom.  The  lowest  two 
(t2u  and  hu)  of  1H  orbitals  are  occupied,  while  the  triply  degenerate 
middle  1H  tiu  orbital  is  the  LUMO  of  superatom.  The  lower  gu 
component  of  2F  orbitials  is  occupied  and  the  higher  t2u  one  is  empty 
as  LUMO  + 1 .  Even  without  closure  of  any  specific  superatom  orbital, 
a  large  HOMO-LUMO  gap  of  1.383  eV  separates  the  occupied  gg 
component  of  the  II  orbital  and  the  empty  tlu  component  of  the  1H 
orbital  of  Sn@Mg12@Sn20.  Therefore,  the  “magic  number”  of  108 
electrons  of  the  icosahedral  matryoshka  clusters  originates  from 
the  high  /#,  symmetry  and  consequently  the  splitting  of  molecular 
orbitals  of  high  angular  momentum  (2F,  1H  and  II)44.  Previously, 
Martin  and  co-workers45  had  shown  that  both  electronic  and  geo¬ 
metric  effects  can  lead  to  extraordinary  stability  for  a  cluster  (i.e.,  a 
magic  cluster). 


Figure  2  |  Energy  levels  (in  eV)  for  the  molecular  orbitals  of  the  Sn@ 
iVlgi2@Sn2()  duster,  which  can  be  assigned  to  a  series  of  superatom 
orbitals  corresponding  to  the  indicated  spherical  harmonics.  A  few 

unoccupied  levels  are  shown  and  highlighted  by  green. 

Meanwhile,  the  suitable  matching  of  atomic  sizes  of  the  compon¬ 
ent  elements  ensures  the  geometry  stability  of  such  unique  three- 
layer  matryoshka  with  high  Iy,  symmetry.  As  summarized  in  the 
Supplementary  Table  S2,  the  covalent  radii46  of  the  constituent  ele- 
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ments  in  an  icosahedral  matryoshka  cluster  have  to  match  with  each 
other  (in  the  range  of  1.22  ~  1.46  A).  Indeed,  we  have  also  investi¬ 
gated  other  possible  icosahedral  matryoshka  A@B12@A2o  clusters 
with  B  =  Be,  Ca  (and  A  =  Sn,  Pb).  However,  all  of  them  exhibited 
several  large  imaginary  vibrational  frequencies  due  to  mismatch  of 
atomic  radius  and  thus  were  not  further  discussed  here. 

The  energy  levels  for  the  four  icosahedral  matryoshka  clusters, 
including  [As@Ni12@As20]3-,  [Sn@Cu12@Sn20]12-,  Sn@Mg12@ 
Sn20,  and  Sn@Zn12@Sn20,  are  compared  in  Supplementary  Figure 
SI.  The  arrangement  of  molecular  orbitals  in  each  species  is  generally 
similar.  The  differences  for  those  clusters  are  related  to  the  specific 
energy  of  each  orbital  as  well  as  the  relative  location  of  the  d  orbitals 
in  the  Ni,  Cu,  or  Zn  derivatives.  However,  the  d  orbitals,  which  locate 
only  in  a  narrow  energy  range,  have  no  noticeable  influence  on  the 
overall  superatom  orbitals.  On  the  other  hand,  the  very  deep  d  levels 
for  As,  Sn,  and  Pb  do  not  overlap  with  the  superatom  orbitals.  As 
discussed  above,  the  universal  arrangement  of  molecular  orbitals 
must  originate  from  the  unique  icosahedral  nesting  doll  structure. 
This  demonstrates  further  the  existence  of  a  class  of  matryoshka 
clusters  with  similar  geometries  and  electronic  properties. 

Analogous  to  Cu,  Ni,  Zn,  Cd  elements  with  filled  d10  shells,  we 
further  design  two  new  icosahedral  matryoshka  clusters  by  using  Mn 
with  a  half-filled  d5  shell,  that  is,  Sn@Mn12@Sn20  and  Pb@Mn12@ 
Pb20.  According  to  Hund’s  rule,  the  Mn  atoms  are  expected  to  carry 
certain  magnetic  moments.  Indeed,  our  DFT  calculations  show  that 
both  Sn@Mn12@Sn20  and  Pb@Mn12@Pb20  possess  a  giant  magnetic 
moment  of  28  pB.  Mulliken  population  analysis  reveals  that  the 
magnetism  mainly  resides  on  the  Mn  atom,  i.e.,  2.880  (3.013)  |rB 
on  each  Mn  atom  for  Sn21Mn12  (Pb2iMn12),  whereas  there  are  cer¬ 
tain  induced  antiferromagnetic  moments  (about  —0.3  ~  0.8  pB)  on 
the  central  and  exterior  Sn  or  Pb  atoms  (see  Table  2).  Even  with  the 
half-filled  d5  shell  and  large  spin  polarization,  these  two  magnetic 
superatoms  exhibit  moderate  HOMO-LUMO  gaps  of  0.382  eV 
(Sn21Mn12)  and  0.614  eV  (Pb21Mn12),  respectively,  which  are  com¬ 
parable  to  the  gap  magnitude  (0.4  ~  0.7  eV)  for  previously  proposed 
magnetic  superatoms,  i.e.,  VNa8-33,  MgsFe34,  Ca8Fe35,  MnCa936, 
MnSr937,  TcMg847,  ScK12  and  ScCsi238. 

The  spin-polarized  density  of  states  for  Sn@Mn12@Sn20  and  Pb@ 
Mn12@Pb20  clusters  are  presented  in  Figure  4.  Even  with  spin  polar¬ 
ization,  the  alignments  of  superatom  orbitals  in  Sn@Mn12@Sn20  and 
Pb@Mn12@Pb20  still  resemble  those  of  their  non-magnetic  counter¬ 
parts  shown  in  Supplementary  Figure  SI.  This  suggests  that  the  3d 
orbitals  of  Mn  in  such  species  are  effectively  non-bonding  and  con¬ 
tribute  only  to  the  overall  magnetic  moment  of  the  clusters.  As 
demonstrated  by  the  spin  density  distribution  in  Figure  4(c)  and 
(d),  the  magnetism  mainly  originates  from  the  unpaired  3d  electrons 
from  the  Mn  atoms,  while  some  superatom  orbitals  (like  2P,  2D,  3S) 
also  exhibit  certain  exchange  splitting  between  the  majority  and 
minority  spin  states.  We  also  observe  some  hybridization  between 
the  3d  states  of  Mn  and  the  superatom  orbitals  (2F,  1G,  1H,  II). 
Previously,  Khanna  and  co-workers32-35  suggested  that  the  hybrid¬ 
ization  between  transition  metal  d  states  and  superatom  states  as  well 
as  the  combined  action  of  the  crystal  field  splitting  and  the  exchange 
splitting  are  the  key  factors  to  stabilize  a  magnetic  superatom  and  to 


Table  2  |  Formation  energy  (Ef),  HOMO-LUMO  gaps  (Ehl),  total 
magnetic  moment  (Mtot)  and  on-site  magnetic  moments  (A,  B,  C 
corresponding  to  Figure  1 )  for  individual  Sn@Mri]  2@Sn2o  and  Pb@ 
Mn  i2@Pb2o  clusters  as  well  as  their  dimers 


Cluster 

Ef  (eV/atom) 

EHL(eV)  M 

ot  ( M-b)  Ma(hb)  Mb  (|.ib) 

Me  (hb) 

Sn?iMm? 

0.346 

0.659 

28 

-0.293  2.880 

-0.698 

Pb2iMni9 

0.507 

0.384 

28 

-0.369  3.013 

-0.781 

[Sn2iMni2]2 

0.341 

0.032 

56 

-0.285  2.869 

-0.716 

[Pb21Mn12]2 

0.503 

0.275 

56 

-0.367  3.009 

-0.779 

ensure  a  large  magnetic  moment.  Although  the  size  and  geometry  of 
the  current  clusters  are  entirely  different  from  the  previous  ones,  this 
design  principle  for  magnetic  superatoms  remains  valid. 

Previously,  an  [Mn@Mn12@Au20]-  cluster,  iso-structural  to  the 
current  icosahedral  matryoshka,  was  also  predicted  to  possess  a  giant 
magnetic  moment  of  44  |iB48,  but  with  a  smaller  HOMO-LUMO  gap 
(0.25  eV).  Using  current  theoretical  scheme,  the  calculated  forma¬ 
tion  energy  for  Mn@Mn12@Au20  is  0.795  eV/atom,  about  twice  of 
Sn21Mn12.  Again,  this  indicates  that  design  of  highly  stable  supera¬ 
toms  requires  the  consideration  of  electronic  structures,  i.e.,  closure 
of  electronic  shell  or  subshell.  Also  note  that  Sun  et  al.49  designed  an 
onion-like  caged  clusters  Fe@Au12@Au42  in  which  the  magnetic 
moment  on  Fe  atom  (3  pB)  is  enhanced  over  bulk  value,  while  Wu 
and  Jena50  recently  proposed  a  series  of  stable  ferromagnetic  hollow 
cages  with  large  magnetic  moments,  i.e.,  Coi2C6  (14  pB),  Mn12C6  (38 
pB),  and  Mn24C18  (70  pB). 

In  addition,  a  [Mn12]  single  molecular  magnet,  i.e., 
Mn12012(CH3C00)i6(H20)4,  which  possesses  a  magnetic  moment 
of  20  pB51,  has  been  intensively  studied52.  To  some  extent,  the  present 
magnetic  superatom  clusters  of  Sn21Mn12  and  Pb21Mn12  can  be 
regarded  as  a  kind  of  ligand-free  single  molecular  magnet  of 
[Mn12]  with  slightly  enhanced  magnetic  moment  (28  pB  vs.  20  pB). 
The  high  symmetries,  appreciable  HOMO-LUMO  gaps,  and  giant 
magnetic  moments  makes  Sn21Mn12  and  Pb21Mn12  clusters  prom¬ 
ising  building  units  for  novel  cluster-assembled  magnetic  materials 
and  devices,  which  have  potential  applications  in  high-density 
information  storage53,  molecular  spintronics54,  quantum  comput¬ 
ing55,  and  magnetic  resonance  imaging  (MRI)56,  etc. 

In  the  cluster-assembled  materials,  clusters  are  either  encapsulated 
into  a  zeolite  matrix57  or  deposited  on  a  substrate58.  It  is  thus  crucial 
to  examine  whether  every  clusters  in  the  assembly  can  keep  their 
identity  without  coalescing  into  larger  clusters,  like  the  previous  case 
ofVCs832.  More  importantly,  the  fantastic  magnetic  characteristics  of 
the  individual  clusters  must  be  well  retained.  As  a  simplest  prototype 
of  cluster  assemblies,  here  we  investigate  the  cluster-cluster  inter¬ 
action  by  considering  the  dimers  of  Sn@Mn12@Sn20  and  Pb@ 
Mn12@Pb20.  The  atomic  structure  and  binding  curve  for  [Pb@ 
Mn12@Pb20]2  dimer  are  shown  in  Figure  5  as  a  representative. 
These  two  clusters  are  stacked  along  a  5-fold  symmetry  axis  with 
parallel  pentagonal  faces  in  a  staggered  orientation,  which  is  more 
stable  than  the  direct  orientation  by  0.180  eV.  Upon  formation  of  the 
dimer,  there  is  no  noticeable  deformation  of  cluster  configuration. 
Accordingly,  the  interaction  strength  between  the  two  spherical  clus¬ 
ters  is  rather  weak,  as  demonstrated  by  the  large  inter-cluster  dis¬ 
tances  (d  =  3.787  A  for  [Sn@Mn12@Sn20]2  and  d  =  4.153  A  for  [Pb@ 
Mn12@Pb20]2,  see  Figure  5)  and  small  inter-cluster  binding  energies 
(0.359  eV  for  [Sn@Mn12@Sn20]2  and  0.291  eV  for  [Pb@Mn12@ 
Pb20]2.  Such  weak  cluster-cluster  interaction  might  be  attributed  to 
the  moderate  HOMO-LUMO  gap  as  well  as  the  sphere-like  ico¬ 
sahedral  geometry  of  the  clusters.  Most  excitingly,  the  total  magnetic 
moments  of  both  cluster  dimers  are  exactly  twice  of  the  individual 
ones,  i.e.,  56  pB-  Indeed,  the  on-site  local  moments  in  the  cluster 
dimers  are  barely  affected  by  the  neighboring  cluster  (see  Table  2). 
To  consider  the  possible  antiferromagnetic  state,  we  initially  set  dif¬ 
ferent  spin  directions  (spin  up  vs.  spin  down)  for  the  Mn  atoms  in 
each  of  the  two  clusters  in  [Sn@Mn12@Sn20]2  or  [Pb@Mn12@Pb20]2 
dimer.  However,  upon  relaxation  of  spin-polarized  electron  wave- 
function,  the  spin  state  of  the  cluster  dimer  spontaneously  trans¬ 
forms  from  antiferromagnetic  to  ferromagnetic.  Moreover, 
inclusion  of  dispersion  correction  for  the  inter- cluster  interaction 
only  moderately  increases  the  magnitude  of  binding  energy  Ein, 
but  has  no  effect  on  the  magnetism  of  the  duster  dimers.  These 
results  clearly  indicate  that  the  Sn@Mn12@Sn20  or  Pb@Mn12@Pb20 
clusters  may  keep  their  identities  and  retain  their  superior  magnet¬ 
ism  in  cluster-assembled  materials  and  devices.  This  is  a  vital  pre¬ 
requisite  for  their  future  applications. 
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Figure  4  j  Spin-polarized  density  of  states  for  (a)  Sn@Mn12@Sn20  and  (c)  Pb@Mn12@Pb20  clusters.  The  occupied  superatom  orbitals  are  labelled. 

The  d  states  from  Mn  atoms  are  highlighted  by  shadow  peaks.  The  HOMO  levels  are  shown  by  dashed  lines.  Atomic  structure  along  with  isosurface  of  spin 
density  distribution  (green  areas)  of  (b)  Sn@Mn12@Sn20  and  (d)  Pb@Mn12@Pb20  clusters  are  also  presented  (color  code:  dark  gray  for  Sn  and  Pb,  light 
purple  for  Mn). 


Conclusion 

A  series  of  three-shell  icosahedral  matryoshka  clusters  related  to  the 
experimentally  known  [As@Ni12@As20]3“  and  [Sn@Cu12@Sn20]12- 
are  proposed.  A  DFT-based  global  minimum  search  confirms  that 
such  icosahedral  matryoshka  structures  are  ground  state  configura¬ 
tions  for  these  A21B12  (A  =  Sn,  Pb;  B  =  Mg,  Zn,  Cd)  binary  clusters. 
In  addition,  the  high  stabilities  of  these  spherical  matryoshka  clusters 
are  demonstrated  by  their  large  HOMO-LUMO  gaps  and  low  forma¬ 
tion  energies.  The  molecular  orbitals  for  different  icosahedral 
matryoshka  clusters  share  the  same  pattern  of  superatom  orbitals. 
The  “magic  number”  of  108  electrons  can  be  attributed  to  the  high  //, 
symmetry  and  consequently  the  splitting  of  superatom  orbitals  of 
high  angular  momentum.  Such  “108-electron  rule”  should  be 
extendable  to  other  icosahedral  matryoshka  superatoms  with  differ¬ 
ent  charge  states,  which  might  lead  to  future  experimental  discov¬ 
eries  of  new  metal  clusters  in  solution  as  well  as  novel  solid-state 
intermetallic  compounds.  Two  related  magnetic  icosahedral 
matryoshka  clusters,  namely  Sn@Mn12@Sn20  and  Pb@Mn12@Pb20, 
are  predicted  to  exhibit  large  magnetic  moments  (28  pB)  and  mod¬ 
erate  HOMO-LUMO  gaps.  More  impressively,  after  formation  of 
cluster  assemblies,  these  two  clusters  are  able  to  keep  their  identities 


and  retain  their  magnetic  moments  due  to  weak  inter-cluster  inter¬ 
action.  These  novel  binary  clusters  with  unique  icosahedral  nesting 
doll  geometry,  high  thermodynamic  stabilities,  and  interesting  phys¬ 
ical  properties  (such  as  optical  gaps  in  the  visible  region  and  giant 
magnetic  moments),  are  expected  to  be  useful  building  blocks  in 
future  nanoscale  materials  and  devices. 

Methods 

Ab  initio  calculations  were  performed  using  the  spin-polarized  density  functional 
theory  as  implemented  in  the  DMol3  program59.  The  core  electrons  were  treated  by  an 
all-electron  relativistic  method  including  scalar  relativistic  effects.  The  generalized 
gradient  approximation  (GGA)  with  PW91  parameterization60  was  adopted  to 
describe  the  exchange-correlation  interaction  and  the  double  numerical  basis  set 
including  d-polarization  functions  (DND)  were  employed.  Vibrational  analyses  were 
performed  for  each  cluster  to  ensure  that  the  optimized  structures  are  the  true  minima 
on  the  potential  energy  surface.  In  addition,  time- dependent  density  functional 
theory  (TD-DFT)  calculations  with  PW91  functional  were  carried  out  for  those 
A@B12@A2o  (A  =  Sn,  Pb;  B  =  Mg,  Zn,  Cd)  clusters  to  obtain  the  optical  adsorption 
spectra  and  optical  gaps. 

To  confirm  that  the  three- shell  icosahedral  configuration  in  Figure  1  is  the  ground 
state  for  the  binary  A21B12  clusters  (A=Sn,  Pb;  B=Mg,  Zn,  Cd,  Mn),  we  performed  an 
unbiased  global  minimum  search  of  two  representative  clusters  (Sn21Mg12  and 
Sn21Zn12)  using  a  genetic  algorithm  (GA)  incorporated  with  DFT  calculations 
implemented  in  the  DMol3  package,  which  was  described  in  our  previous  publica- 
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Figure  5  |  Atomic  structure  (upper)  and  binding  curve  (lower)  for  [Pb@ 
Mn12@Pb2ol2  dimer.  The  equilibrium  distance  (d)  and  optimal  inter¬ 
cluster  binding  energy  (Eint)  are  labelled. 


tion61.  In  the  GA  search,  sixteen  initial  configurations  were  generated  from  scratch. 
Any  two  individuals  in  this  population  were  then  chosen  as  parents  to  produce  a  child 
cluster  via  a  “cut  and  splice”  crossover  operation,  followed  by  an  optional  mutation 
operation  (either  a  small  random  displacement  on  each  atom  or  exchange  of  a  pair  of 
different  types  of  atoms)  of  35%  probability.  The  diversity  of  the  populations  was 
filtered  by  the  inertia  of  each  cluster.  After  about  3000  GA  iterations,  the  icosahedral 
matryoshka  structure  was  obtained  for  both  Sn@Mg12@Sn20  and  Sn@Zn12@Sn20 
clusters  as  the  ground  state.  In  addition,  some  metastable  isomers  (Supplementary 
Figure  S2)  are  found  from  GA-DFT  search,  but  are  less  energetically  favorable  and 
will  not  be  further  discussed. 
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